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Intermittency transition to generalized synchronization in coupled time-delay systems 
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In this paper, we report the nature of transition to generaUzed synchronization (GS) in a system 
of two coupled scalar piecewise linear time-delay systems using the auxiliary system approach. We 
demonstrate that the transition to GS occurs via on-off intermittency route and also it exhibits 
characteristically distinct behaviors for different coupling configurations. In particular, the inter- 
mittency transition occurs in a rather broad range of coupling strength for error feedback coupling 
configuration and in a narrow range of coupling strength for direct feedback coupling configuration. 
It is also shown that the intermittent dynamics displays periodic bursts of period equal to the de- 
lay time of the response system in the former case, while they occur in random time intervals of 
finite duration in the latter case. The robustness of these transitions with system parameters and 
delay times has also been studied for both linear and nonlinear coupling configurations. The results 
are corroborated analytically by suitable stability conditions for asymptotically stable synchronized 
states and numerically by the probability of synchronization and by the transition of sufcLyapunov 
exponents of the coupled time-delay systems. We have also indicated the reason behind these dis- 
tinct transitions by referring to unstable periodic orbit theory of intermittency synchronization in 
low-dimensional systems. 

PACS numbers: 05.45.Xt,05.45.Pq 



I. INTRODUCTION 

Synchronization of interacting chaotic oscillators is one 
of the most interesting nonlinear phenomenon and is an 
inherent part of many natural systems (cf . P, Q ) • The 
concept of synchronization is receiving a central impor- 
tance in recent research in nonlinear dynamics due to its 
potential applications in diverse areas of science and tech- 
nology- Since the identification of chaotic synchroniza- 
tion ^ 0, Q , several works have appeared in identifying 
and demonstrating basic kinds of synchronization both 
theoretically and experimentally [1, 2, 3, 4, 5]. There are 
also attempts to find a unifying framework for defining 
the overall class of chaotic synchronizations [1, 0j Hi • 

One of the interesting synchronization behaviors of 
unidirectionally coupled chaotic systems is the general- 
ized synchronization (GS), which was conceptually in- 
troduced in Ref. Q. Generalized synchronization is ob- 
served in coupled nonidentical systems, where there ex- 
ists some functional relation between the drive X{t) and 
the response Y{t) systems, that is, Y{t) = F{X{t)). 
With GS, all the response systems coupled to the drive 
lose their intrinsic chaoticity (sensitivity to initial condi- 
tions) under the same driving and follow the same tra- 
jectory Hence the presence of GS can be detected us- 
ing the so called auxiliary system approach 10], where 
an additional system (auxiliary system) identical to the 
response system is coupled to the drive in similar fash- 



ion. Auxiliary system approach is particularly appealing 
since it can be implemented directly in an experiment 
and, in addition, this method allows one to utilize ana- 
lytical approaches for studying GS. However, one has to 
be aware that if there are multiple basins of attraction 
for the coupled drive-response system, then the auxiliary 
system approach can fail. 

Generalized synchronization (GS) has been well stud- 
ied and understood in systems with few degrees of free- 
dom and for discrete maps [1, M, [O, d, El, Q [ll- 
The concept of GS has also been extended to spatially ex- 
tended chaotic systems such as coupled Ginzburg-Landau 
equations [l6|. Recently, the terminology intermittent 
generalized synchronization (IGS) ^l^] was introduced 
in diffusively coupled Rossler systems in analogy with 
intermittent lag synchronization (ILS) 
termittent phase synchronization (IPS) 
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18 , 11911 and in- 
2a,lM[2i, and 

also experimentally in coupled Chua's circuit. Very re- 
cently, it has been shown [23| that the transition to inter- 
mittent chaotic synchronization (in the case of complete 
synchronization) is characteristically distinct for geomet- 
rically different chaotic attractors. In particular, it has 
been shown that for phase coherent chaotic attractors 
(Rossler attractor) the transition occurs immediately as 
soon as the coupling strength is increased from zero and 
for non- phase-coherent attractors (Lorenz attractor), the 
transition occurs slowly as the coupling is increased from 
zero. 

Time-delay systems form an important class of dynam- 
ical systems and recently they are receiving central im- 
portance in investigating the phenomenon of chaotic syn- 
chronizations in view of their infinite dimensional nature 
and feasibility of experimental realization [2^. [25l. [26l. [27j . 
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FIG. 1: The hyperchaotic attractor of the system ([T]) for the 
parameter values a — 1.0, fo = 1.2 and r = 15.0. 



While the concept of GS has been well established in 
low dimensional systems, it has not yet been studied in 
detail in coupled time-delay systems and only very few 
recent studies have been dealt with GS in time-delay sys- 
tems [23, [25||. In particular, the mechanism of onset of 
GS in coupled time-delay systems and its characteristic 
properties have not yet been clearly understood and re- 
quire urgent attention. 

In this paper, we investigate the characteristic proper- 
ties of nature of onset of GS from asynchronous state in 
unidirectionally coupled piece- wise linear time-delay sys- 
tems exhibiting highly non-phase coherent hyperchaotic 
attractors [2^] . We find that the onset of GS is preceded 
by on-off intermittency mechanism from the desynchro- 
nized state. We have also identified that the intermit- 
tency transition to GS exhibits characteristically distinct 
behaviors for different coupling schemes. In particular, 
the intermittency transition occurs in a broad range of 
coupling strength for error feedback coupling configura- 
tions and in a narrow range of coupling strength for direct 
feedback coupling configurations, beyond certain thresh- 
old value of coupling strength. In addition, the intermit- 
tent dynamics is characterized by periodic bursts away 
from the temporal synchronized state with period equal 
to the delay time of the response system in the case of 
broad range intermittency transition whereas it is char- 
acterized by random time intervals in the case of narrow 
range intermittency transition. We have also confirmed 
these dynamical behaviours in both linear and nonlinear 
coupling configurations. We have analyzed these transi- 
tions analytically using Krasvoskii-Lyapunov functional 
approach and numerically by the probability of synchro- 
nization and by the sw&Lyapunov exponents. We have 
also addressed the reason behind these transitions using 
periodic orbit theory. The robustness of these transi- 
tions with the system parameters in both the linear and 
nonlinear, error feedback and direct feedback coupling 
configurations, are also studied. 

The plan of the paper is as follows. In Sec. II, we 
will point out the existence of broad range intermittency 
transition to GS for linear error feedback coupling pro- 
portional to {xi{t) — X2{t)) while in Sec. Ill the exis- 
tence of narrow range intermittency transition is shown 
for the linear direct feedback coupling of the form xi (t) , 
where xi{t) and X2{t) are the drive and response sig- 
nals, respectively (see below for details). In Sec. IV we 



will discuss the existence of broad range intermittency 
transition for nonlinear error feedback coupling of the 
form (/(.Ti(i — T2)) — f {x2{t — T2))) , where f[x) is an odd 
piece- wise linear function. The existence of narrow range 
intermittency transition is discussed in Sec. V for nonlin- 
ear direct feedback coupling of the form {f{xi{t — T2)). 
Finally in Sec. VI, we will summarize our results. 



II. BROAD RANGE (SLOW/DELAYED) 
INTERMITTENCY TRANSITION TO GS FOR 
LINEAR ERROR FEEDBACK COUPLING OF 
THE FORM {xi{t) - X2{t)) 

We consider the following first order delay differential 
equation introduced by Lu and He ^28] and discussed in 
detail by Thangavel et al. (29j . 



x{t) = -ax{t) +bf{x{t -t)), 



(1) 



where a and b are parameters, t is the time-delay and / 
is an odd piecewise linear function defined as 



0, a: < -4/3 

-1.5a; -2, -4/3 < a; < -0. 

f{x) = { X, -0.8 < a; < 0.8 

-1.5X + 2, 0.8 < a; < 4/3 

0, a; > 4/3. 



(2) 



Recently, we have reported [13] that the system ([!]) 
exhibits hyperchaotic behavior for suitable parametric 
values. For our present study, we find that for the choice 
of the parameters a = 1.0,6 = 1.2 and r = 15.0 with 
the initial condition x{t) = 0.9, t e (—5,0), Eq. ^ ex- 
hibits hyperchaos. The corresponding pseudoattractor is 
shown in Fig. [TJ The hyperchaotic nature of Eq. ^ is 
confirmed by the existence of multiple positive Lyapunov 
exponents. The first ten maximal Lyapunov exponents 
for the above choice of parameters as a function of delay 
time r is shown in Fig. [5] (spectrum of Lyapunov expo- 
nents in this paper are calculated using the procedure 
suggested by Farmer fsTl). 

To be specific, we first consider the following unidi- 
rectional, linearly coupled systems with drive xi{t), re- 
sponse X2{t) and an auxiliary x^lt), 

xi{t) = - axi{t) + bif{xi{t - Ti)), (3a) 
X2{t) = - ax2{t) + b2f{x2(t - T2)) 

+ b3{xi{t)-X2{t)), (3b) 

xait) = - axsit) + b2f{x3{t - T2)) 

+ b3{xi{t)-X3{t)), (3c) 

where 61, 62 and 63 are constant parameters, and ri and 
T2 are constant delay parameters. Note that when bi ^ 62 
or Ti T2 or both, corresponding to parameter mis- 
matches, we have unidirectionally coupled nonidentical 
systems (Eq. ([3^) and dSb)), while the auxiliary system 
is given by ([31;) and f{x) is the odd piece- wise linear 
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FIG. 2: The first ten maximal Lyapunov exponents \max of 
tfie scalar time-delay system for the parameter values a = 
1.0, h\ = 1.2, r e (2, 29) (which is the same as Eq. HJ with 6i 
replaced by 6) ((3^). 



function The coupling in ^p) may be also called a 
linear error feedback coupling. 

For simplicity, we have chosen h\ = 62 so that the 
time-delays t\ and T2 alone introduce a simple form of 
parameter mismatch between the drive x\ (t) and the re- 
sponse X2{t). We have chosen the values of parameters 
as a = 1.0,61 = 62 = 1-2, n = 20 and T2 = 25. For 
this parametric choice, in the absence of coupling, all 
the three systems ^ evolve independently and exhibit 
hyperchaotic attractors, which is confirmed by the exis- 
tence of multiple positive Lyapunov exponents (Fig. [2]). 
The actual value of the positive Lyapunov exponents 
for r = 20 are 0.00916, 0.00759, 0.00565, 0.00283 and 
0.00073 and for r 25 they are 0.01234, 0.01067, 
0.00886, 0.00658, 0.00386, 0.00229, 0.00123 and 0.00033. 



A. Stability Condition 

With GS, as all the response systems under the same 
driving follow the same trajectory, it is sufficient to iden- 
tify the existence condition for establishment of complete 
synchronization (CS) between the response X2{t) and the 
auxiliary X3{t) systems in order to achieve GS between 
the drive Xi{t) and the response X2{t) systems. 

Now, for CS to occur between the response X2{t) and 
the auxiliary X3{t) variables, we consider the time evo- 
lution of the difference system with the state variable 
A — X3{t) — X2{t). It can be written for small values of 
A as 



where 



-{a + b3)A + b2f{x2{t-T2))Ar,, (4) 



-1.5, -4/3 <x < -0.8 
f'{x) = { 1, -0.8 < a; < 0.8 (5) 
-1.5, 0.8<x<4/3. 



FIG. 3: Dynamics in the phase space of the systems (|3]). (a) 
and (b): Approximate GS and CS, respectively, for the value 
of the coupling strength &3 = 0.4. (c) and (d): Perfect GS 
and CS, respectively, for the value of the coupling strength 
63 = 0.9. 



The synchronization manifold, X2(t) = x^it), is locally 
attracting if the origin of this equation is stable. Fol- 
lowing Krasovskii-Lyapunov functional approach, we de- 
fine a positive definite Lyapunov functional of the form 
[27l . [33 . i33)] (details of stability analysis are given in ap- 
pendix El 



J To 



(6) 



where /i is an arbitrary positive parameter, /i > 0. 
The solution of Eq. namely A = 0, is stable if 
the derivative of the functional along the trajectory of 
Eq. (HI) is negative. This negativity condition is satisfied 
if 63 -|- a > ''^^ i^2it-r2)) _|_ from which it turns out 
that the sufficient condition for asymptotic stability is 



a + bs> \b2f(x2{t-T2))\. 



(7) 



Now from the form of the piecewise linear function f{x) 
given by Eq. we have. 



\f{^2{t~T2))\ = 



1.5, 0.8 < \x2\ < I 
1.0, |x2|<0.8. 



(8) 



Consequently the stability condition ([7]) becomes a-f 63 > 
1 1.562 1 > 1 62 1- Thus one can take 



a + 63 > I62 



(9) 



as the less stringent or approximate stability condition 
(as the synchronization dynamics of the coupled systems 
can occur even beyond the inner region, \x2\ < 0.8) 
for ^ to be valid, while 



a-l-63 > 11.562 



(10) 
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FIG. 4: The intermittent dynamics of the response X2{t) and 
auxiliary xz{t) systems for the value of the coupling strength 
&3 = 0.4. (a) Time traces of the difference X2{t) — X2.(t) corre- 
sponding to Fig.[3}3, (b) Enlarged in x scale to show bursts at 
periodic intervals when bursts of larger amplitudes A > |0.01| 
are considered and (c) Enlarged in y scale to show random 
bursts when bursts of smaller amplitudes A < |0.01| are con- 
sidered. 



can be considered as the most general or stringent or 
exact stability condition (as the full synchronization dy- 
namics of the coupled systems ^ lies within the outer 
region 0.8 < \x2\ < |) specified by ([7]) for asymptotic 
stability of the synchronized state A = 0. 



B. Approximate (Intermittent) Generalized 
Synchronization 

In order to understand the mechanism of transition 
to synchronized state, it will be important to follow the 
dynamics from the parameter values at which the less 
stringent condition is satisfied. Figure [3^ shows the ap- 
proximate GS (which may also be termed as intermit- 



0.2 
0.15 
$ 0.1 



0.05 



(a) 



^ 




< 



0.1 



0.01 



0.001 



le-04 



20 40 60 80 100 
t 

(b) ; 




0.1 



10 



FIG. 5: The statistical distribution of laminar phases corre- 
sponding to the Fig.|4l (a) For threshold value A = |0.1| and 
(b) for the threshold value of A = |0.0001|. 



tent generalized synchronization (IGS) in analogy with 
intermittent lag synchronization (ILS)) between the drive 
xi{t) and the response X2{t) systems, whereas Fig. [5Jd 
shows the approximate CS between the response X2(t) 
and the auxiliary x^ (t) systems for the values of the pa- 
rameters a = 1.0, 6i = 62 = 1.2, Ti = 20, T2 — 25 and 
63 = 0.4 satisfying the less stringent condition Per- 
fect GS and perfect CS are shown in Figs. [3j: and [3]i 
respectively for 63 = 0.9 satisfying the general stability 
condition (fTU]) . Time traces of the difference X2 (t) — X3 (t) 
corresponding to approximate CS (Fig. are shown in 
Fig. m which show periodic bursts with period between 
two consecutive bursts approximately equal to the time- 
delay of the response system t « 25 when 'on' states 
of amplitude greater than |0.01| are considered. Fig. 2)3 
shows an enlarged (in x scale) part of Fig. [4^ to view 
the bursts at periodic intervals when bursts of larger am- 
plitudes (A > |0.01|) are considered, while Fig. Uj; is an 
enlarged (in y scale) version of Fig. [4)d to show random 
bursts when bursts of smaller amplitude A < |0.01| are 
considered. 

Usually the intermittent dynamics is characterized by 
the entrainment of the dynamical variables in random 
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time intervals of finite duration [s^, [s^. But from 
Fig. it is evident that the intermittent dynamics dis- 
plays periodic bursts from the synchronous state with 
period approximately equal to the delay time of the re- 
sponse system, when amplitudes of the state variable 
A| — \x3{t) — X2it)\ > 0.01 are considered, for the val- 
ues of the coupling strength at which the less stringent 
stability condition ^ is satisfied. The statistical fea- 
tures associated with the intermittent dynamics is ana- 
lyzed by calculating the distribution of laminar phases 
A{t) with amplitude less than a threshold value to A. A 
universal asymptotic power law distribution A{t) oc t~" 
is observed for the threshold value A — 1 0.0001 1 with the 
value of the exponent a = —1.5 as shown in Fig. [SId, 
which is quite typical for on-off intermittency. On the 
other hand the distribution of laminar phases A{t) for 
the value of the threshold value of |A — 0.1| shows a 
periodic structures (Fig. [5^), where the peaks occur ap- 
proximately ai t = nT, n — 1,2,..., where T « T2 is the 
period of the lowest periodic orbit of the uncoupled sys- 
tem ([Hb)- Note that —3/2 power law is observed for the 
intermittent dynamics shown in Fig.[3]for laminar phases 
A{t) with amplitude less than A = |0.0001| (as an illus- 
trative example), which is also evident from Fig.|4|:, while 
periodic bursts are observed for 'on' state of amplitude 
greater than |0.01|. It is to be noted that such periodic 
bursts of period approximately equal to the time-delay 
of the response system for larger threshold value of A 
along with the on-off intermittency behavior for larger 
threshold value of A has also been observed by Zhan et 
al in unidirectionally coupled Mackey- Glass systems 
, where the authors discussed relation between two modes 
of synchronization, namely, CS and GS. 



C. Characterization of IGS 

Now we characterize the intermittency transition to GS 
by using (i) the notion of the probability of synchroniza- 
tion $(63) as a function of the coupling strength &3 [23| . 
which can be defined as the fraction of time during which 
\x2(t) — x^{t)\ < e occurs, where e is a small but arbi- 
trary threshold, and (ii) from the changes in the sign 
of sM&Lyapunov exponents (which are nothing but the 
Lyapunov exponents of the subsystem) in the spectrum 
of Lyapunov exponents of the coupled time-delay sys- 
tems. Fig. shows the probability of synchronization 
$(63) as a function of the coupling strength 63 calcu- 
lated from the variables of the response X2{t) and the 
auxiliary 2:3 (i) systems for CS between them. For the 
range of 63 G (0,0.39), there is absence of any entrain- 
ment between the systems resulting in asynchronous be- 
havior and the probability of synchronization $(63) is 
practically zero in this region. However, starting from 
the value of 63 = 0.39 and above, there appear oscil- 
lations in the value of the probability of synchronization 
$(63) between zero and some finite values less than unity, 
exhibiting intermittency transition to GS in the range of 
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FIG. 6: (a) The probability of synchronization ^{bs) between 
the response X2{t) and the auxiliary 2:3 (t) systems and (b) 
Largest Lyapunov exponents of the coupled drive xi{t) and 
response X2{t) systems and (|3Jd). 



63 € (0.4, 0.62) for which the less stringent stability con- 
dition ([9]) is satisfied. Beyond 63 = 0.62, $(63) attains 
unit value indicating perfect GS. Note that the above in- 
termittency transition occurs in a rather wide range of 
the coupling strength (this can also be termed as slow 
or delayed intermittency transition in analogy with the 
terminology used in 23]), which has also been confirmed 
from the transition of successive largest sit6Lyapunov 
exponents in the corresponding range of the coupling 
strength. 

The spectrum of the first fifteen largest Lyapunov ex- 
ponents Xmax of the coupled drive xi (t) and response 
2:2 (t) systems are shown in Fig. ^jp. From the general 
stability condition it is evident that for the cho- 
sen value of the parameter a = 1.0, the less stringent 
stability condition Q is satisfied for the values of cou- 
pling strength 63 > 0.2. Correspondingly, the least posi- 
tive sw6Lyapunov exponent of the response system ^p) 
gradually becomes negative from 63 > 0.2. Subsequently, 
the remaining positive sw&Lyapunov exponents gradually 
become negative and attain saturation in the range of 
63 e (0.2,0.8). This is in accordance with the fact that 
the less stringent stability condition is satisfied only in 
the corresponding range of coupling strength 63. This 
is a strong indication of the broad range intermittency 
(IGS) transition to GS. For 63 > 0.8, the general sta- 
bility condition (jTU]) is satisfied, where one can observe 
perfect GS as is evidenced by both the probability of syn- 
chronization approaching unit value and by the negative 




FIG. 7: The probability of synchronization $(63) in 3- 
dimensional plots (a) as a function of the system parameter 62 
and the coupling strength 63 and (b) as a function of the cou- 
pling delay T2 and the coupling strength 63, exhibiting broad 
range intermittency transition to GS for linear error feedback 
coupling. 



saturation of sitfeLyapunov exponents calculated between 
the drive and response systems. The inference is that the 
correlation between the oscillations of the systems even- 
tually becomes stronger with the strength of the cou- 
pling, and this is indicated by the successive transition 
of sw6Lyapunov exponents to negative values. 

It is a well established fact that a chaotic attractor 
can be considered as a pool of infinitely many unsta- 
ble periodic orbits of all periods. Synchronization be- 
tween two coupled systems is said to be asymptotically 
stable, if all the unstable periodic orbits of the response 
system are stabilized in the transverse direction of the 
synchroniation manifold. Consequently, all the trajecto- 
ries transverse to the synchronization manifold converge 
for suitable values of coupling strength and this is re- 
flected in the negative values of the transverse Lyapunov 
exponents (sM&Lyapunov exponents) upon synchroniza- 
tion (23j . From our results, we find that the sw^Lyapunov 
exponents gradually become negative in a broad range of 
coupling strength 63 after certain threshold value and 
this is in accordance with the known results on gradual 
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FIG. 8: Dynamics in the phase space of the systems (|lip : (a) 
and (b) Approximate GS and CS, respectively, for the value 
of the coupling strength 63 — 0.64 and (c), (d) Perfect GS 
and CS, respectively, for the value of the coupling strength 
h = 0.8. 



stabilization of unstable periodic orbits of the response 
system in the complex synchronization manifold of low 
dimensional systems p3| . Unfortunately, methods for lo- 
cating UPO's and calculating their transverse Lyapunov 
exponents have not been well established for time-delay 
systems and hence a quantitative proof for the gradual 
stabilization of UPO's has not been given here. How- 
ever, the gradual stabilization of UPO's along with their 
transverse Lyapunov exponent in the range of intermit- 
tency transition have been reported for the case of cou- 
pled Rossler and Lorenz systems in Ref. [1^ . It can then 
be inferred from these studies that the broad range inter- 
mitteny transition in the case of error feedback coupling 
configuration is due to the fact that the strength of the 
coupling 63 contributes only less significantly to stabilize 
the UPO's as the error xi{t) — X2{t) gradually becomes 
smaller from the transition regime after certain threshold 
value of the coupling strength. 



The robustness of the intermittency transition in a 
broad range of coupling strength with the system pa- 
rameter 62 G (1.1,1.6) and with the coupling delay 
T2 e (10, 20) has also been confirmed. Fig. [7^ shows the 
3-dimensional plot of the probability of synchronization 
as a function of the system parameter 62 and the coupling 
strength 63, while Fig. [Tb shows the 3-dimensional plot 
of $(63) as a function of the coupling delay T2 and the 
coupling strength 63. The above figures (Fig. [7]) clearly 
reveal the broad range intermittency transition to GS in 
the case of linear error feedback coupling scheme. 
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III. NARROW RANGE (IMMEDIATE) 
INTERMITTENCY TRANSITION TO GS FOR 
LINEAR DIRECT FEEDBACK COUPLING OF 
THE FORM xi{t) 

To illustrate the narrow range intermittency transition 
to GS, we consider the unidirectional linear direct feed- 
back coupling of the form 

±i{t) = -axi{t) + bJ{xi{t-T,)), (11a) 
i2it) = -aa;2(t) + fe2/(x2(t-T2)) + 63x1(0, (lib) 

X3{t) = -aX3{t) + b2f{x3{t-T2)) + b3Xi{t),{llc) 

where f{x) is of the same odd piece- wise linear form as in 
Eq. ([2]) . Assuming the same values of the parameters as 
before and proceeding in the same way as in the previous 
case, one can obtain the sufficient condition for asymp- 
totically stable CS between the response X2{t) (lllb ) and 
the auxiliary 0:3 (t) pTb ) systems as 

a>\b2f'{x2it-T2))\. (12) 

It is to be noted that the above stability condition 
holds good only for the case when coupling is present, 
that is &3 ^ 0. When there is no coupling (63 = 0), by 
definition, there will be a desynchronized chaotic state. 
As soon as the value of the coupling strength is increased 
from zero, the stability condition always lead to syn- 
chronized state even for very feeble values of 63 for pa- 
rameters satisfying the stability condition, as it is inde- 
pendent of the coupling strength 63. As the values of the 
parameters satisfying the stability condition (|12p rapidly 
leads to immediate transition to synchronized state as 
soon as the coupling is switched on, it is difficult to iden- 
tify the possible transitions to synchronized state. In 
addition, as the stability condition is independent of the 
coupling strength 63, one is not able to explore the dy- 
namical transitions as a function of coupling strength for 
the parameter values satisfying the stability condition 
(|12p . Hence we study the synchronization transition by 
choosing the parameters violating the stability condition 
as a = 1.0,61 = 1.2,62 = 1-1, Ti = r2 = 20 and varying 
the coupling strength 63 in order to identify the mech- 
anism of synchronization transition. Here, in this case 
61 and 62 alone introduce the parameter mismatch while 
Ti ~ T2 (It may be added that the qualitative nature of 
the dynamical transitions remain the same even when the 
mismatch is either in time delays alone, that is ti ^ T2, 
61 = 62 or in both the system parameters and time de- 
lays, 61 7^ 62, Ti 7^ T2, as confirmed below in the three 
dimensional plots of Figs. [T^ . 

As 63 is varied from zero, transition from desynchro- 
nized state to approximate GS occurs for 63 > 0.6. Ap- 
proximate GS (IGS) between the drive system xi{t) spec- 
ified by Eq. pTb ). and the response system X2{t) repre- 
sented by Eq. pTb ). is shown in Fig. [8^ whereas the 
approximate CS between the response X2{t) (Eq. dHh)) 
and the auxiliary X3{t) (Eq. (|llb )) systems is shown in 
Fig. [SJd for the value of the coupling strength 63 = 0.64. 
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FIG. 9: The intermittent dynamics of the response X2{t) (|llb ) 
and auxiliary (t) <\llt ) systems for the value of the coupling 
strength 63 = 0.64. (a) and (b) Time traces of the difference 
X2{t) — X3{t) corresponding to Fig. [SId. 



Prefect GS and CS are shown in Figs.[8j: and[8ji, respec- 
tively, for the value of the coupling strength 63 = 0.8. 
The intermittent dynamics at the transition regime cor- 
responding to the value of the coupling strength 63 — 0.64 
is shown in Figs. [51 in which Fig. shows the enlarged 
part of Fig. [5^. It is clear from this figure that the inter- 
mittent dynamics displays intermittent bursts at random 
time intervals. The statistical distribution of the laminar 
phases again shows a universal asymptotic —1.5 power 
law behavior for the threshold value A = |0.0001|, which 
is typical for on-off intermittency transitions, as shown 
in Fig. [lOl 

Now we characterize the intermittency transition to 
GS in the present case, again by using the notion of 
the probability of synchronization $(63) and from the 
changes in the sign of sM6Lyapunov exponents of the 
coupled system. The probability of synchronization is 
shown in Fig. [Tib as a function of the coupling strength, 
again calculated from the response 2:2 (t) and the auxil- 
iary X3{t) systems, Eqs. (|llb ) and ([Tib ), respectively, 
which remains zero in the range of 63 G (0, 0.60) and 
oscillates between its maximum and minimum values in 
a narrow range of 63 £ (0.60,0.68) confirming the ex- 
istence of approximate CS in the later range. Above 
63 — 0.68 the probability of synchronization acquires 
its maximum value depicting perfect CS between the 
response X2(t) and the auxiliary X3(t) systems. Cor- 
respondingly there exists perfect GS between the drive 
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FIG. 10: The statistical distribution ol laminar phases for the 
Fig. [9] 



xi{t) and the response X2{t) systems. Fig. [TTb shows 
the first twelve maximal Lyapunov exponents of the cou- 
pled drive xiit) and the response x^it) systems. The 
least positive sM^Lyapunov exponent of the response sys- 
tem starts to become negative from 63 > 0.60. Subse- 
quently, all the other positive sw&Lyapunov exponents 
become negative and reach saturation in a rather nar- 
row range of 63 £ (0.60,0.68). Thus the narrow range 
intermittency (IGS) transition (this can also be termed 
as immediate intermittency transition in analogy with 
the terminology used in j23| ) is confirmed from both 
the probability of synchronization, calculated from the 
response and the auxiliary systems, and negative satura- 
tion of SM&Lyapunov exponents, calculated from the drive 
and the response systems. 

As discussed in the previous section, the narrow range 
intermittency transition is in accordance with the stabi- 
lization of all the unstable periodic orbits of the response 
system in a narrow range as a function of the coupling 
strength and this is reflected in the immediate transi- 
tion of all the SM^Lyapunov exponents (Fig. [TTb) to neg- 
ative values. It is also to be noted that the narrow range 
intermitteny transition in the case of direct feedback cou- 
pling configuration can be attributed to the fact that the 
strength of the coupling 63 contributes significantly pro- 
portional to the strength of the signal xi{t) to stabilize 
all the UPO's immediately at the transition regime after 
certain threshold value of the coupling strength as in the 
case of low-dimensional systems [23!| . 

The robustness of the intermittency transition in a 
narrow range of the coupling strength 63 for a range 
of values of the parameter 62 G (1.1,1.6) and the delay 
Ti = T2 G (10, 20) is shown in Fig. [T^ The 3-dimensional 
plot of the probability of synchronization as a function of 
the system parameter 62 and the coupling strength 63 is 
shown in Fig.fT^. while Fig. 112b shows the 3-dimensional 
plot of $(63) as a function of the coupling delay T2 and 
the coupling strength 63. 
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FIG. 11: (a) The probability of synchronization $(63) be- 
tween the response X2{t) (|llb) and the auxiliary x^it) Ullb l 
systems and (b) Largest Lyapunov exponents of the coupled 
drive x\(t) and response X2{t) systems (|llb ,) and (|llb V 



IV. BROAD RANGE INTERMITTENCY 
TRANSITION TO GS FOR NONLINEAR ERROR 
FEEDBACK COUPLING OF THE FORM 

(f[x^{t-r2))- f{x2{t-r2))) 

Next we demonstrate the existence of the above types 
of distinct characteristic transitions for nonlinear cou- 
pling configurations as well. For this purpose, we con- 
sider the unidirectional nonlinear error feedback coupling 
of the form 

xi{t) = - axiit) + bif{xi{t - Ti)), (13a) 
X2{t) = - ax2it) + b2j{x2{t - T2)) 

+ h^{f{xi{t - T2)) - f{x2{t - r2))), (13b) 
X3.{t) = - aa;3(t) + b2f{x'i{t - T2)) 

+ h{f{xi{t - T2)) - f{x^{t - T2))), (13c) 

where f{x) is again of the same piece- wise linear form as 
in Eq. The parameters are now fixed as a = 1.0, 5i = 
62 = 1.2, Ti — 20 and the coupling delay T2 — 25, where 
the delays alone form the parameter mismatch between 
the drive and response systems in Eqs. (|13p . Following 
Krasvoskii-Lyapunov theory, for complete synchroniza- 
tion so that the manifold A = ^3 {t) — X2 (t) between the 
response X2{t) and the auxiliary x^{t) approaches zero, 
one can obtain the stability condition as 



a>\{b2-b:,)f{x2{t^T2))\. 



(14) 
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FIG. 12: The probability of synchronization ^(fes) in 3- 
dimensional plots (a) as a function of the system parameter 

62 and the coupling strength 63 and (b) as a function of the 
coupling delay T2 and the coupling strength 63, exhibiting 
narrow range intermittency transition to GS for linear direct 
feedback coupling 

Consequently from the form of the piecewise linear 
function the stability condition (|14p becomes a > 
11.5(62 — 63)1 > 1(62 — 63)1. Thus one can take 

a> \b2- hi (15) 

as less stringent condition for ([7]) to be valid, while 

a > 11.562 - l-Sfoal (16) 

can be considered as the most general condition speci- 
fied by (|14p for asymptotic stability of the synchronized 
state A = X2{t) — 2^3 (t) = 0. For the chosen values 
of the parameters, the less stringent stability condition 
is satisfied for the values of the coupling strength 

63 e (0.2, 0.535) and the general stability condition (fT6|) 
is satisfied for 63 > 0.535. 

As the coupling strength is increased from zero, ap- 
proximate GS occurs from 63 > 0.2. Figure [T5b shows 
the approximate GS (IGS) between the drive xi{t) 
(Eq. p^ll and the response X2{t) (Eq. (fT^ )) systems 
for the value of the coupling strength 63 = 0.37, while the 
approximate CS between the response X2{t) (Eq. ([I3h)) 
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FIG. 13: Dynamics in the phase space of the systems (I13|l 'l. 
(a) and (b): Approximate GS and CS, respectively, for the 
value of the coupling strength 63 = 0.37. (c) and (d): Perfect 
GS and CS, respectively, for the value of the coupling strength 
63 = 0.8. 



and the auxiliary X3{t) (Eq. (|13b)) systems is shown in 
Fig. (HSb). Perfect GS and perfect CS are shown in 
Figs.[T5b and[T2Ji respectively for 63 = 0.8. The intermit- 
tent dynamics exhibited by the coupled systems at the 
transition regime is shown in Fig. 1141 which shows bursts 
at the period approximately equal to the delay time of 
the response system X2 (t) for bursts of amplitude greater 
than |0.01| (Fig.fHb). The statistical distribution of the 
laminar phases away from the intermittent bursts shows 
an asymptotic —1.5 power law behavior for the thresh- 
old value A = |0.0001| (see Fig. [H;), typical for on-off 
intermittency, which is shown in Fig. 115b . On the other 
hand for the threshold value A = |0.1| Fig.fTSb shows pe- 
riodic structures similar to Fig. [S^ with peaks occuring 
approximately at t = nT,n = 1,2, where T « r2 is 
the period of the lowest periodic orbit of the uncoupled 
system ([T3b)- 

Now, the intermittency transition is again character- 
ized using the probability of synchronization and the 
sM^Lyapunov exponents as in the previous cases. Fig. 116b . 
shows the probability of synchronization $(63), the value 
of which remains zero in the range 63 G (0, 0.2) due to 
the fact that there lacks any entrainment between the re- 
sponse X2(t) and the auxiliary X3{t) systems, whereas it 
fluctuates between the two extreme values in a rather 
broad range of the coupling strength 63 e (0.2,0.42), 
depicting the existence of intermittency transition in 
the corresponding range of 63. Perfect CS exists for 
63 > 0.42 as evidenced from the maximum value of 
<i>(63). Correspondingly there exists perfect GS between 
the drive Xi{t) and the response X2{t) systems. Figure 
116b shows the transition of sufoLyapunov exponents of 
the spectrum of Lyapunov exponents of the coupled drive 
xi{t) (Eq. (113^)) and the response X2{t) (Eq. (|13b)) sys- 
tems. The su6Lyapunov exponents become negative in 
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FIG. 15: The statistical distribution of laminar phases cor- 
responding to the Fig. 1141 (a) For threshold value A = |Q.1| 
and (b) for the threshold value of A = |0.0001|. 



FIG. 14: The intermittent dynamics of the response 12 (t) 
(|13b) and auxiliary X3{t) (|13b ) systems for the the value of 
the coupling strength 63 = 0.37 for nonlinear error feedback 
coupling, (a) Time traces of the difference 2:2 (t) — 2:3 (t) corre- 
sponding to Fig. 113b . (b) Enlarged in x scale to show bursts at 
periodic intervals when bursts of larger amplitudes A > |0.01 
are considered and (c) Enlarged in y scale to show random 
bursts when bursts of smaller amplitudes A < |0.Q1| are con- 
sidered. 



the range 63 £ (0.2, 0.42) confirming the broad range in- 
termittency (IGS) transition in a rather wide range of 
the couphng strength and this is again due to the grad- 
ual stabilization of the unstable periodic orbits of the 
response systems because of the less significant contri- 
bution of the coupling strength 63 as the error becomes 
gradually smaller from the transition regime beyond cer- 
tain threshold value of the coupling strength as discussed 
in Sec. II. The robustness of the intermittency transi- 
tion with the system parameter &2 and the coupling de- 
lay T2 as a function of coupling strength 63 is shown as 
3-dimensional plots in Figs. [T71 



V. NARROW RANGE INTERMITTENCY 

TRANSITION TO GS FOR NONLINEAR 
DIRECT FEEDBACK COUPLING OF THE 
FORM f{xi{t-T2)) 

Now we consider the unidirectional nonlinear coupling 
of the form 

xi{t) = - axi{t) + bif{xi{t ~ n)), (17a) 

X2{t) = - ax2{t) -f b2f{x2{t - T2)) 

+ b3f{x,it-T2)), (17b) 
X3{t) = - ax3{t) + b2f(xz{t - T2)) 

+ b3f{xi{t-T2)). (17c) 

Choosing the values of the parameters as in the pre- 
vious case and following Krasvoskii-Lyapunov functional 
approach for the asymptotically stable synchronized state 
A = X3{t) — X2it) — 0, one can obtain the sufficient condi- 
tion for asymptotic stability for complete synchronization 
of the response X2 (t) and the auxiliary x^ (t) systems as 

a>\b2f'{x2{t-T2))\. (18) 



11 



1.2 
1 

0.8 
& 0.6 
0.4 
0.2 




0.02 



(a) 



0.01 



-0.01 



-0.02 



T 



0.2 



0.4 



0.6 



0.8 




FIG. 16: (a) The probability of synchronization $(63) be- 
tween the response X2{t) (|13b) and the auxiliary X3{t) <\13t ) 
systems and (b) Largest Lyapunov exponents of the coupled 
drive xi{t) and response X2{t) systems (|13b ,) and (|13b V 
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The above stability condition rapidy leads to immedi- 
ate transition to synchronized state even for very feeble 
values of the coupling strength 63 for the parameter val- 
ues satisfying the stability condition as the stability 
condition is independent of 63 as in the previous linear 
coupling case (Sec. IIIip . Hence it is difficult to identify 
the possible dynamical transitions to synchronized state 
as a function of the coupling strength 63. Hence we study 
the synchronization transition as a function of the cou- 
pling strength 63 by choosing the parameters violating 
the stability condition as a = 1.0, 61 = 1.2, &2 = 1-1 and 

Tl = T2 = 15. 

As 63 is varied from zero for the above values of the 
parameters, transition from desynchronized state to ap- 
proximate GS occurs for 63 > 0.74. The approximate 
GS (IGS) between the drive xi{t) and the response X2{t) 
variables described by Eqs. lTTk ) and (|17b ) is shown in 
Fig. [TSk . whereas Fig. [TSb shows the approximate GS 
between the response X2(t) and the auxiliary X3{t) vari- 
ables (Eqs. ^TTh ) and (fT7b )) for the value of the cou- 
pling strength 63 — 0.78. Perfect GS and perfect GS are 
shown in Figs. [TSb and [TSH respectively for 63 = 0.9. 
Time traces of the difference X2 (t) — X3 {t) corresponding 
to approximate GS (Fig. [TSb) is shown in Fig. [121 which 
shows intermittent dynamics with the entrainment of the 
dynamical variables in random time intervals of finite 
duration. Fig. 119b shows the enlarged picture of part 
of Fig. [TOb . The statistical distribution of the laminar 
phases again shows a universal asymptotic —1.5 power 



FIG. 17: The probability of synchronization $(63) in 3- 
dimensional plots (a) as a function of the system parameter 
62 and the couphng strength 63 and (b) as a function of the 
coupling delay T2 and the couphng strength 63 for the case 
of nonlinear error feedback coupling configuration given by 
Eq. (|13|) . showing broad range intermittency transition. 
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FIG. 18: Dynamics in the phase space of the systems (I17|l l. 
(a) and (b): Approximate GS and CS, respectively, for the 
value of the coupling strength 63 — 0.78. (c) and (d): Perfect 
GS and CS, respectively, for the value of the coupling strength 
63 = 0.9. 
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FIG. 19: The intermittent dynamics of the response X2{t) 
(|17b ) and auxihary 3:3 (t) (117b ) systems for the the value of 
the couphng strength 63 — 0.78. (a) and (b) Time traces of 
the difference 2:2 (t) — 2:3 (t) corresponding to Fig. 118b . 



FIG. 21: (a) The probabihty of synchronization $(63) be- 
tween the response X2{t) (|17b) and the auxihary 3:3 (t) ((T7J;) 
systems and (b) Largest Lyapunov exponents of the coupled 
drive xi{t) and response X2{t) systems (|17b ,) and (|17b V 
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FIG. 20: The statistical distribution of laminar phase for the 
Fig. [191 



law behavior for the threshold value A — |0.0001|, typi- 
cal for on-off intermittency, as shown in Fig. [501 

As in the previous cases, now we characterize the in- 
termittency transition to GS using the notion of prob- 
ability of synchronization ^(b^) and from the changes 
in the sign of sw^Lyapunov exponents in the spectrum 
of Lyapunov exponents of the coupled time-delay sys- 
tems ((T7)) . Fig. \Tlh shows the probability of synchro- 
nization $(63) as a function of the coupling strength 



63 calculated from the response X2{t) and the auxiliary 
xsit) variables (Eqs. p7b ) and (HZt)) for CS between 
them. In the range of 63 G (0,0.74), the probability of 
synchronization remains approximately zero. Upon in- 
creasing the value of 63, $(^3) oscillates in the narrow 
range of 63 G (0.74, 0.78) depicting the existence of in- 
termittency transition. This narrow range transition is 
also confirmed from the transition of successive largest 
sw6Lyapunov exponents. The spectrum of the first nine 
largest Lyapunov exponents Xmax of the coupled drive 
xi{t) and response X2{t) variables (Eqs. (|17b .) and (|17b )) 
is shown in Fig. [2Tb. It is also evident from the spectrum 
that the su&Lyapunov exponents suddenly become neg- 
ative in the narrow range of 63 G (0.74,0.78), and then 
reach saturation values for 63 > 0.78. This confirms the 
narrow range intermittency (IGS) transition to GS. This 
is also in accordance with the immediate stabilization 
of all the UPO's of the response system as discussed in 
Sec. ([ml) . 

The robustness of the transition with the system pa- 
rameter 62 and the delay time T2 as a function of coupling 
strength 63 is shown as 3-dimensional plots in Fig. [211 

VI. CONCLUSION 

In conclusion, we have studied intermittency transition 
to generalized synchronization from desynchronized state 
in unidirectionally coupled scalar piecewise linear time- 
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are analyzed using the analytical stability condition for 
synchronized state, probability of synchronization $(63) 
between the response and the auxiliary systems, and by 
the changes in the sign of sw&Lyapunov exponents in the 
spectrum of Lyapunov exponents of the drive and the 
response systems in both the linear and nonlinearly cou- 
pled time-delay systems. In spite of the fact that both 
the probability of synchronization and the sw^Lyapunov 
exponents have been calculated from different systems, 
we have found good agreement between them in showing 
intermittency transition in all the cases. We hope this 
study will contribute to the basic understanding of the 
nature of transition to GS in coupled time-delay systems 
and we are now investigating the experimental verifica- 
tion of these findings in nonlinear electronic circuits. 
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APPENDIX A: STABILITY CONDITION 



FIG. 22: The probability of synchronization ^(fes) in 3- 
dimensional plots (a) as a function of the system parameter 
62 and the coupling strength 63 and (b) as a function of the 
coupling delay T2 and the coupling strength 63 for the case 
of nonlinear direct feedback coupling configuration given by 
Eq. (I17|) showing narrow range intermittency transition to 
GS. 

delay systems for different coupling configurations using 
the auxiliary system approach. We have shown that the 
intermittency transition to GS occurs in a broad range 
of coupling strength for both linear and nonlinear error 
feedback coupling configurations whereas it occurs in a 
narrow range of coupling strength for both the linear and 
nonlinear direct feedback coupling configuration. It has 
also been shown that the intermittent dynamics displays 
periodic intermittent bursts of period equal to the delay 
time of the response systems in the former cases and it 
takes place in random time intervals in the latter cases. 
The robustness of the intermittent dynamics with the 
system parameters and the delay time is also studied as 
a function of the coupling strength for both error feed- 
back and direct feedback (linear and nonlinear) coupling 
confi gur ations. Universality of these intermittent behav- 
iors [36| (periodic and random) and their (broad and 
narrow range) transitions are also confirmed for differ- 
ent forms of linear and nonlinear coupling configurations 
in other time-delay systems such as Mackey-Glass and 
Ikeda systems. These distinct characteristic behaviors 



To estimate a sufficient condition for the stability of the 
solution A = 0, we require the derivative of the functional 
V{t) along the trajectory of Eq. ([4]), 

^ = -(a + &3)A2 + 62/'(.T2(t - T2))AA,, + - ^A^^, 

(Al) 

to be negative. The above equation can be rewritten as 
^ = -/iA2r(X,M), (A2) 

where X = A,JA, T =[{{a + bs - fj.)/ fj.)-{b2f'ix2it - 
T2))/y)X + X'^]. In order to show that ^ < for 
all A and At-j and so for all X, it is sufficient to show 
that Tmin > 0. One can easily check that the absolute 
minimum of F occurs at X — ■^b2f'{x2{t — T2)) with 

r™i„ =[4/i(a + 63 -/i) -blf'{x2{t-T2)f]/A^?. Conse- 
quently, we have the condition for stability as 

^, , blr{x2{t-T2)f , ^, , 

a + bz> hM = *(M)- (A3) 

4^ 

Again ^(/i) as a function of ji for a given f'{x) has an ab- 
solute minimum at /i = (|62/'(a^2(^— T'2))|)/2 with $mm = 
\h2f{x2(t - T2))|. Since $ > = |&2/'(x2(t - T2))|, 

from the inequality (jA3|) . it turns out that the sufficient 
condition for asymptotic stability is 

a + b^>\b2f\x2{t-T2))\. (A4) 
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Now from the form of the piecewise Hnear function f{x) 
given by Eq. [21 we have, 

ill <ol^ 

Note that the region |x2| > 4/3 is outside the dynamics 
of the present system (see Eq. ([2])). Consequently the 



stabihty condition (|A4p becomes a + 63 > 1.5|62| > \b2\. 

Thus one can take a + 63 > I&2I as a less stringent 
condition for (|A4p to be valid, while 

a>1.5|&2|, (A6) 

as the most general condition specified by (|A4[) for 
asymptotic stability of the synchronized state A = 0. 
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